Strong enhancement of the valley splitting in a 2D electron system in silicon 
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Using magnetocapacitance data, we directly determine the chemical potential jump in a strongly 
correlated 2D electron system in silicon when the filling factor traverses the valley gap at v — 1 and 
v = 3. The data yield a valley gap that is strongly enhanced compared to the single-particle value 
and increases linearly with magnetic field. This result has not been explained by existing theories. 
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A two-dimensional (2D) electron gas in (lOO)-silicon 
metal-oxide-semiconductor field-effect transistors (MOS- 
FETs) is a unique double-layer electron system with 
strong interlayer interactions. Indeed, the valley index 
is identical with the isospin quantum number, and one 
can therefore try to apply to this system a good deal 
of recent theoretical work on double-layer electron sys- 
tems. In a double-layer system with small layer sepa- 
ration compared to the interelectron distance in cither 
layer, strong interlayer correlations are predicted to give 
rise to the appearance of both novel ground states and 
novel excitations [jl], ||. In quantizing magnetic fields 
strong enough to fully polarize real electron spins, two 
cases as determined by the isospin system symmetry are 
discussed: (i) at finite layer separation the symmetry is 
U(l) corresponding to an "easy-plane" anisotropy and 
(ii) at vanishing layer separation the symmetry trans- 
forms to SU(2) and the anisotropy disappears. In the 
first (second) case the lowest energy charge-carrying ex- 
citations are merons (skyrmions) that are isospin textures 
with the charge e/2 (e), where e is the electron charge 
jl], g]. There is little doubt that the valley splitting in 
the 2D electron system in silicon MOSFETs should be 
of many-body origin at least for lowest odd filling factors 
because the electron-electron interactions (and correla- 
tions) in this system are strong; particularly, in accessi- 
ble magnetic fields, the Coulomb energy exceeds signifi- 
cantly the cyclotron energy. However, this is the strongly 
interacting limit in which existing theories are not valid 
and, therefore, they cannot be directly applied to silicon 
MOSFETs. The origin of the excitations for the valley 
splitting is unknown so far. 

Experimental investigations of the valley splitting were 
performed largely (§) at high filling factor v > 9 based 
on analysis of the beating pattern of Shubnikov-de Haas 
oscillations in tilted magnetic fields |J, H H 01 ■ The gap 
value, its linear dependence on substrate bias and 
its insensitivity to parallel magnetic field j^] are consis- 
tent with single-particle theoretical considerations for an 
asymmetric potential well that contains a 2D electron gas 
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In this paper, we perform, for the first time, low- 
temperature measurements of the chemical potential 
jump across the valley gap at the lowest filling factors 
v = 1 and v = 3 in a 2D electron system in silicon using 
a magnetocapacitance technique. The valley splitting is 



found to exceed strongly the single-particle value, decay- 
ing with filling factor. Unexpectedly, the data are best 
described by a linear increase of the valley gap with mag- 
netic field, which is similar to the proportional magnetic 
field dependence of the enhanced spin gap for the 2D 
electrons in AlGaAs/GaAs heterostructures jij], O]. 

Measurements were made in an Oxford dilution refrig- 
erator with a base temperature of w 30 mK on high- 
mobility (lOO)-silicon MOSFETs (with a peak mobility 
close to 2 m 2 /Vs at 4.2 K) having the Corbino geometry 
with diameters 250 and 660 /jm. The gate voltage was 
modulated with a small ac voltage 15 mV at frequencies 
in the range 2.5 - 25 Hz and the imaginary current com- 
ponent was measured with high precision using a current- 
voltage converter and a lock-in amplifier. Care was taken 
to reach the low frequency limit where the magnetoca- 
pacitance is not distorted by lateral transport effects. A 
dip in the magnetocapacitance at integer filling factor is 
directly related to a jump of the chemical potential across 
a corresponding gap in the spectrum of the 2D electron 
system [ 14], |]| . 

Typical traces of the magnetocapacitance, C(B), at 
different electron densities, n s , are displayed in Fig. [|. 
It oscillates as a function of filling factor, v — hcn s /eB, 
reflecting the modulation of the thermodynamic density 
of states (DOS) in quantizing magnetic fields. Narrow 
dips in the magnetocapacitance at integer v are separated 
by broad maxima where the value C approaches in the 
high field limit the geometric capacitance, Co, between 
the gate and the 2D electrons, in agreement with previous 
studies El: 
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where A is the sample area, D — dn s /dfj, is the thermo- 
dynamic DOS, and fi is the chemical potential. We have 
checked that in the range of electron densities used, the 
value D(B = 0), which is calculated from the magnetoca- 
pacitance in the high field limit (Co) and that in B = 0, 
corresponds to 2m/irH 2 , where m = 0.19m e and m e is 
the free electron mass. Note that at the high field edge 
of the dip in C(B), the capacitance C can overshoot Co 
forming a local maximum (so-called negative compress- 
ibility effect the effect being not so pronounced at 
v = 3 as at v — 1, see Fig. [j]. 
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FIG. 1: Magnetocapacitance in the low frequency limit for 
electron densities 5.3 x 10 11 cm~ 2 (a) and 2.2 x 10 11 cm -2 
(b). The level of the geometric capacitance Co is indicated by 
a dotted line. 



Integrating Eq. ([!]) over a dip at integer v — i>q yields 
a jump, A, of the chemical potential between two neigh- 
boring quantum levels J15[ 
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where the integration over B is equivalent to the one 
over n s as long as the dip is narrow. We note that, ex- 
perimentally, it is easier to analyze C(B) traces: while 
being independent of magnetic field, the geometric capac- 
itance Co increases slightly with n s as the 2D electrons 
are forced closer to the interface. 

Apparently, if the capacitance C does not reach Co 
(i.e., the maximum thermodynamic DOS is insufficiently 
large), the jump A will be smaller than the level split- 
ting by the level width. While the jump A is still de- 
termined by the area of the dip in C(B) (see the in- 
set of Fig. ||(a)), we estimate the level splitting includ- 
ing the level width contribution with the help of a sim- 
ilar integration of Eq. (|l|) between the magnetic fields 
B\ = hcn s /e(vo + 1/2) and B2 = hcn s /e(vo — 1/2). 

In Fig. |^(a), we show the v = 3 minimum in the C(B) 
curve for n s = 3.85 x 10 11 cm~ 2 at different tempera- 
tures. As the temperature, T, is lowered, the minimum 
deepens appreciably and becomes narrower until it satu- 
rates in the low T limit, which is consistent with temper- 
ature smearing of the thermodynamic DOS. In a similar 
way, the maximum capacitance decreases with temper- 
ature departing from the geometric capacitance Co- In 
insufficiently high magnetic fields, the (positive) differ- 
ence between Co and C at the low field edge of the dip 
is larger compared to that at the high field edge. There- 
fore, we determine the jump A by replacing the reference 
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FIG. 2: The temperature dependence of the v — 3 minimum 
in the magnetocapacitance at n s = 3.85 x 10 11 cm~ 2 (a) and 
of the chemical potential jump with (squares) and without 
(circles) the level width contribution (b). The dashed lines 
in (b) are guides to the eye. The inset shows the same mini- 
mum in C(B) along with the shifted reference level (solid and 
dashed lines) used to determine the jump A of Eq. (^) when 
the geometric capacitance Co (dotted line) is not attained. 



level Co in Eq. (g) by the capacitance value at the low 
field edge of the dip as depicted in the inset to Fig. ||(a). 
This gives a little underestimated values of A for such 
minima whose weakly asymmetric shape persists more 
or less, irrespective of magnetic field. 

The chemical potential jump A(T) across the v = 3 
valley gap for the data of Fig. ||(a) with and without al- 
lowing for the temperature-dependent level width is dis- 
played in Fig. |2|(b) . The value A increases with decreas- 
ing T and saturates in the low temperature limit, this 
temperature dependence being weaker compared to that 
of the minimum magnetocapacitance (and the minimum 
D). The estimated valley splitting including the level 
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FIG. 3: The valley gap for v = 1 (circles), v = 3 (squares), 
and v = 5 (triangles) as a function of magnetic field (a) and 
electron density (b). The systematic error bars correspond to 
the level width contribution. A linear fit of the v — 3 data is 
shown by a dashed line. The solid line in (b) is the expected 
linear n 3 dependence of the single-particle valley gap |j| for 
depletion layer charge density 1 x 10 11 cm -2 . 

width contribution is practically independent of temper- 
ature, as expected. It is clear that the experimental un- 
certainty of determination of the valley splitting is the 
smallest in the low temperature limit where the modula- 
tion of the thermodynamic DOS is maximal. 

In Fig. we show the v = 1 and v = 3 valley splitting 
measured in the low temperature limit as a function of 
both magnetic field and electron density. For both filling 
factors, the value of valley gap turns out to be strongly 
enhanced compared to the single-particle value |J (solid 
line in Fig. |](b)). The valley gap enhancement decays 
rapidly with v so that almost no enhancement is observed 
for v = 5. To our surprise, the data are best described by 
a linear increase of the valley gap with B (or n s ). This 



fact is reliably established for v — 3 where the geometric 
capacitance Co is attained in the entire range of magnetic 
fields used except for the lowest B, which is indicated 
by systematic error bars in Fig. [| that correspond to 
the level width contribution. Extrapolation to B = of 
the linear magnetic field dependence of the v = 3 gap 
yields a value which is consistent with the single-particle 
splitting. For filling factor v = 1, the dependence of A 
on B (or n s ) can also be described by a linear function, 
although the experimental uncertainty (estimated from 
data dispersion) of the data for v = 1 is markedly larger 
than that for v = 3. 

The linear magnetic field dependence of the enhanced 
valley gap observed in silicon MOSFETs is similar to 
the puzzling proportional B dependence of the enhanced 
spin gap in the 2D electron system in AlGaAs / GaAs het- 
erostructures jl2], Before making comparison with 
existing theories, we would like to emphasize that they 
are expected to hold in the limit of weak electron-electron 
interactions, i.e., the Coulomb energy, E c — e 2 /el (where 
I = (hc/eB) 1 / 2 is the magnetic length), in a 2D electron 
system should be small compared to the cyclotron energy, 
huj c . For real 2D electron systems, one determines that in 
a magnetic field of 10 T, the Coulomb and the cyclotron 
energies are approximately equal to each other for the 2D 
electrons in GaAs, whereas in Si MOSFETs the energy 
E c exceeds hjj c by a factor of 4. Therefore, even for the 
case of GaAs, the validity of the Landau-level-based con- 
siderations jy, ^, O, |l^] is questionable. Formally, these 
cannot be applied to the case of Si MOSFETs. Bearing 
this in mind, we will make qualitative comparison of our 
results with the aforementioned theories below. 

Theoretically, the basic idea is that the gap en- 
hancement is expected to be controlled by electron- 
electron interactions whose strength is characterized by 
the Coulomb energy E c [jl], g, [ll], 0. Hence, as long as 
the enhancement dominates the single-particle gap value, 
Ao, the many-body enhanced gap should essentially be 
proportional to B 1 / 2 . This law is in contradiction to the 
experiment. Not to mention that the predicted values 
of the many-body gap are an order of magnitude larger 
than the experimental ones. 

Two approaches to introduce corrections to the B 1 ! 2 
law have been formulated. Based on the traditional 
theory of exchange-enhanced gaps JD] ], the authors of 
Ref. (17) took account of the corrections to the exchange 
energy due to level overlap and finite 2D layer thickness 
as well as the correlation energy contribution. Knowing 
that in our case the level overlap is small [jl3|, we have 
verified that the approach |[7| is unable to noticeably in- 
crease the power of the theoretical B x l 2 dependence of 
the gap. 

In accordance with the other approach, in the limit 
E c 3> Ao, a possible formation of the isospin textures 
- skyrmions or merons 0] - - is expected to lead 
to a reduction of the exchange-enhanced gap which is 
related to single isospin-flip excitations. We have verified 
that in the range of magnetic fields (or electron densities) 
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used, the predicted crossover between the two regimes as 
governed by the ratio Aq / E c [[[) has a weak effect on the 
theoretical square-root dependence A(-B). Although the 
actual symmetry of the isospin system in Si MOSFETs 
is unknown, it is unlikely that the outcomes for the two 
textures Jl], ||] would be very different, because the layer 
separation is small. 

Thus, for our case both approaches fail to modify ap- 
preciably the B 1 / 2 law and, therefore, they are not even 
able to give a qualitative account of our experimental 
data. 

In summary, we have performed low-temperature mea- 
surements of the chemical potential jump across the val- 
ley gap at the lowest v ~ 1 and v = 3 in a 2D electron 
system in silicon. The valley splitting is found to be 
strongly enhanced compared to the single-particle value, 
decaying with filling factor. The data are best described 
by a linear increase of the valley gap with magnetic field. 
This result is similar to the proportional B dependence 



of the enhanced spin gap for the 2D electrons in GaAs 

Recently, in the limit E c hiu c , a remarkably differ- 
ent behavior of the many-body gap to create a charge- 
carrying spin texture excitation at integer v has been 
predicted jl8|: A st = v\Q\hu> c , where Q is the integer 
topological charge. If the exchange effects play a similar 
role for the valley splitting |9| , the predicted linear B de- 
pendence of the gap is consistent with our findings. Still, 
both the gap value and its dependence on filling factor 
are not explained by the theory fl8|| . 
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